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Abstract 

We compute, in topological terms, the spectral flow of an arbitrary family of self-adjoint 
Dirac type operators with classical (local) boundary conditions on a compact Riemannian 
manifold with boundary under the assumption that the initial and terminal operators of the 
family are conjugate by a bundle automorphism. This result is used to study conditions for 
the existence of nonzero spectral flow of a family of self-adjoint Dirac type operators with 
local boundary conditions in a two-dimensional domain with nontrivial topology. Possible 
physical realizations of nonzero spectral flow are discussed. 

Keywords: Aharonov-Bohm effect, massless Dirac fermions, graphene, topological insu- 
lators, self-adjoint Dirac operator, spectral flow, Atiyah-Singer index theorem, Atiyah-Bott 
index theorem, index locality principle. 

1 Introduction 

Not to mention high-energy physics and quantum field theory, the ideas of modern geometry and 
topology become increasingly important in condensed matter physics [THE]- In particular, the 
Atiyah-Singer index theorem [10] explains a topological protection of zero-energy Landau level 
and related peculiarities of the quantum Hall effect in graphene [QE]. Topologically protected 
zero modes play an essential role in the motion of vortices in superfluid helium-3 [SUIT]. The 
quantum Hall effect 0[3] and topological insulators [8,9| are examples of the states of matter 
with topological order parameter. 

The Aharonov-Bohm effect (ABE) [12,13] has actually initiated this development. A mag- 
netic flux localized in a region completely unavailable for a quantum particle (e.g., surrounded 
by infinitely high potential barrier) nevertheless affects its motion, modifying the geometry of 
quantum space. A periodic dependence of electron energy levels in a ring as a function of the 
magnetic flux through the ring resulting in appearance of a persistent current (see Ref. [13] and 
references therein) is a bright manifestation of ABE. When the change of the magnetic flow is 
equal to an integer number of the flux quanta, the energetic spectrum should return to its initial 
state. Until recently, ABE was studied mainly for usual nonrelativistic electrons described by 
the Schrodinger equation. After discovery of graphene, the ABE for ultrarelativistic electrons 
described by Dirac equation with zero mass has attracted attention |15H17j . From the mathe- 
matical point of view, this is a much richer case. The Dirac operator is not semibounded and 
hence its spectral flow [18] can be nonzero. It is worth noting that the nonzero spectral flow 
of the Dirac operator has been discussed already in a context of condensed matter physics. It 
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Figure 1: Example of a domain X with m = 5 boundary components. The bold lines show the 
boundary components d + X = dXiUdX^UdX^ on which B > in the boundary condition in ([T|). 
(It is only these components that contribute to the spectral flow according to formula ([3]).) The 
arrows show the positive sense of going around the contour when computing the winding number 
of the gauge transformation \i along the corresponding boundary component. 

results in additional forces ("Kopnin forces") acting on vortices in superfluid helium-3 [5jlllj. 
Coming back to ABE, it means that the coincidence of the whole energy spectra at the change of 
the magnetic flux at integer number of the flux quanta does not necessarily mean the periodicity 
of individual eigenenergies (like the shift m — > m + 1 transforms Z to itself). Nonzero spectral 
flow corresponds to a physical situation when an adiabatically slowly varying magnetic field leads 
to a production of electron-hole (or, in general, particle-antiparticle) pairs from the vacuum: 
"positron" levels cross zero-energy level transforming into electron ones. The vacuum reconstruc- 
tion effects were discussed in physics of superfluid helium-3 [5] and in physics of graphene [7] 
but without any relation with ABE. Here we study conditions of existence of nonzero spectral 
flow of a family of Dirac-like self-adjoint operators with local boundary conditions in a domain 
with nontrivial topology. 

In a bounded domain IcM 2 with smooth boundary dX (see Fig. [1]), consider the boundary 
value problem 



where n x and n y are the inward normal components and B is a nonvanishing real- valued function 
on the boundary. (Berry and Mondragon [19] were the first to consider boundary conditions 
of this kind.) The operator Dq corresponding to this problem is self-adjoint and Fredholm on 
L 2 (X, C 2 ). Next, let /i be a smooth function on X with = 1. The "gauge transformation" 
Dq — > [iDq^ 1 takes Dq to the operator D\ = Dq + Qi (where Q± is a self-adjoint matrix 
function) with the same boundary conditions. In physicists' language, the gradient of the phase 
[i is an abelian (U(l)) gauge field, i.e., an electromagnetic vector potential. We consider a two- 
dimensional domain with m — 1 holes pierced by magnetic flux tubes. The case in which all 
magnetic fluxes through the holes are integer multiples of the magnetic flux quantum corresponds 
to [i = 1. Let us join Q\ with Qq = by a continuous family Q t , t £ [0, 1], of self-adjoint matrix 
functions. The spectral flow sfjL^} of the family 




in X, (n y — in x )ui = Bui on dX, (1) 



A = D Q + Q t 



(2) 
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i.e., the number of eigenvalues of Dt that changed their sign from minus to plus as the parameter 
t varies from to 1 minus the number of eigenvalues that changed their sign from plus to minus, 
does not change under continuous deformations of the family provided that Dq and D\ remain 
isospectral in the course of deformation. As far as the authors know, the problem of finding 
this spectral flow (also for families of Dirac operators of more general form) was posed for the 
first time and partially solved in [20], where the spectral flow was computed up to an integer 
factor c m depending on the number m of boundary components. Further, it was shown in [20] 
that C2 = 1, and it was conjectured that c m = 1 for all m. In the present paper, we establish 
a general result (see Theorem [2] below), which, in particular, proves this conjecture to be true. 
Thus, Theorem 1 in [20j acquires the following form: 

Theorem 1. The spectral flow of the family ([2]) is given by the formula 



is the winding number of the restriction of the function \jl to d + X. ( The set d + X is a union 
of finitely many circles; when defining the winding number, the positive sense of any of these 
circles is the one for which the domain X remains to the left when moving along the circle.) 

This theorem shows that the coefficients c m that remained unfound in |20|, Theorem 1] are 
equal to unity for all m. The same is true for [2Q|, Theorems 2 and 3]; all unknown coefficients 
c m occurring there are equal to unity. 

Theorem [1] follows from a general result established in the present paper. We give a computa- 
tion in topological terms of the spectral flow of an arbitrary family {Dt}, t £ [0, 1], of self-adjoint 
Dirac type operators with local boundary conditions on a compact Riemannian manifold X with 
boundary dX under the assumption that D\ = UDqII^ 1 , where U is some automorphism of 
the bundle in which the operators Dt act. (In contrast with [20J, we assume neither that the 
principal part of Dt is independent of t nor even that the principal parts of Do an d D\ coincide.) 
Namely, we prove (see Theorem [2] below) that 



where the right-hand side is the index of an elliptic operator with boundary conditions on the 
manifold X x S 1 with boundary, t being the coordinate on the circle S 1 and the operator D 
being obtained from the family {Dt} by clutching the operators Dq and D\ with the use of 
the automorphism U. (Recall that formula (HJ) for families of self-adjoint elliptic operators on 
a closed manifold X was established in [IB]-) The right-hand side of ^ can be computed by 
the Atiyah-Bott formula [2]] (see also [22J Sec. 20.3]). Note, however, that we do not rely 
on the Atiyah-Bott formula in the proof of Theorem HJ relation Q between the spectral flow 
and the index permits one to use the localization method (see [23Tf25] ) and cut the domain 
into parts, thus reducing the problem to the case of a domain with one hole (m = 2), for 
which a straightforward computation was carried out in [20] . Note also that the localization 
method proves to be an important technical tool when proving relation Q itself. The proof is in 
many aspects similar to that in [261 Proposition 5.6] of the spectral flow formula for families of 
differential operators Agranovich-Vishik elliptic with parameter on a closed compact manifold 
but contains a number of new important lines of argument related to the presence of boundary 
conditions. 



sf D t = wind u, 

d+x 



(3) 



where d + X is the part of dX where B > and 





(4) 
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Figure 2: Definition of spectral flow. There are three eigenvalues \j(t), j = 1,2,3, contributing 
to the spectral flow. Both the computation according to the definition and counting the number 
of zero crossings (with regard to direction) give the value —1 for the spectral flow. 

2 Spectral Flow 

Recall the definition of spectral flow in the form presented in [26J (cf. |27j). Let {Bt}, t G [0, 1], 
be a family, continuous in the sense of uniform resolvent convergence, of unbounded self-adjoint 
operators with purely discrete spectrum on a Hilbert space T~L. Then there exists a partition 
= to < t\ < t2 < • • • < t n+ i = 1 of the interval [0, 1] and real numbers 71, ... , 7 n +i such that 
jj does not lie in the spectrum Spec(-Bi) of the operator B t for t £ [tj-\,tj], 71 = 7 n +i < 0, and 
if 71 < 0, then the half-open interval [71,0) does not contain any points of spectrum of Bq and 
Bi. 

Definition 1 (see [26] . Definition A. 18). The spectral flow of the family {Bt}, t £ [0, 1], is the 
numbe 

n 

sf{B t } = ^ m i si g n (7j - Ij+i), (5) 

where rrij is the number of eigenvalues (counting multiplicities) of the operator Bt j on the 
interval between 7^ and 7j+i- 

This definition is illustrated in Fig. [21 which, in particular, clarifies why this definition is 
consistent with the notion of spectral flow as the number of eigenvalues passing through zero 
(with direction taken into account). 

3 Main Results 

Let E be an even-dimensional Hermitian vector bundle over a compact Riemannian manifold X 
with boundary, and let 

A : C°° (X, E) — ► C°° (X, E) (6) 

1 - ) The right-hand side of formula ^ is independent of the choice of the partition {tj} and the numbers jj by 
Theorem A. 19 in [26]. 
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be a formally self-adjoint Dirac type operatoro Next, let a subbundle L C Ey of dimension 
dim L = I dim E be given in the restriction Ey of the bundle E to the boundary Y = dX of 
the manifold X such that 

((TA^nCx))^) ±L X VxGY, (7) 

where L x is the fiber of L at x and n(x) is the unit inward conormal vector on the boundary. 
Consider the operator ([6]) on the set of sections u £ C°° (X, E) satisfying the homogeneous 
boundary condition 

^l(u\y) = 0, where ir^: Ey — > Ey/L is the natural projection. (8) 

(In other words, u(x) £ L x for x £ Y.) In particular, the boundary condition in ([T]) is of 
this form. It is well known (see [29] and [28} Chaps. 18 and 19]) that the boundary condition 

is elliptic, the operator ^ with domain given by this condition is essentially self- adjoint 
on L 2 (X,E), and its closure Al is an unbounded Fredholm self-adjoint operator on L 2 (X,E) 
with discrete spectrum and with domain consisting of sections u belonging to the Sobolev space 
H 1 (X,E) and satisfying condition (JSj) (in which u\y is treated as the element of H l / 2 (Y,Ey) 
obtained from u by restriction to Y by virtue of the trace theorem and ttl is treated as a mapping 
tt l : HV\Y,E Y ) — > H l l 2 {Y,E Y /L)). 

Now assume that both the Dirac type operator A ([6]) and the subbundle L continuously 
depend on a parameter t G [0, 1] (namely, the coefficients of A and ttl depend on t continuously 
together with all of their derivatives^]); i.e., A = A(t) and L = L(t). Moreover, assume that 
condition ([7]) holds for each t. Then, by Theorem 7.16 in [SO], the operator A(t) L ^ continuously 
depends on t in the topology of uniform resolvent convergence, and Definition [1] specifies the 
spectral flow si{A(t) L u\} of the family {A(t) L ^}, t G [0, 1]. 

Next, let an automorphism U : E — > E of the bundle E be given such that 

A(l) = UA(0)U-\ U(L{0)) = L(l). (9) 

Then U A{Q) L ^U~ l = A(l) L ^y, i.e., the operators A(0)u \ and ^(l)^) are similar and hence 
isospectral, so that the spectral flow of the family {A(t)u t \} is a homotopy invariant (in the 
class of families satisfying a condition of the form ©). Thus, it is natural to pose the problem 
of computing it in topological terms. 

To do this, we introduce an auxiliary elliptic boundary value problem on the Cartesian 
product X x 5 1 of the manifold X by the circle S 1 (see Fig. [3|). Namely, let us define a bundle 
£ over X x S 1 as follows. Take the pullback of E to the product X x [0, 1] via the natural 
projection 7r : X X [0, 1] — >• X and then use the automorphism U : (it* E) Xx {o} — ^ ( 7r *^)xx{i} 
as the clutching automorphism Q By conditions Q, the family {^4(i)} specifies a well-defined 
differential operator on the space of sections of the bundle £, while the family of subbundles L(t) 
defines a subbundle C C Syxs 1 m the restriction of £ to the boundary Y x S 1 of the manifold 
X x S 1 . 

Proposition 1. The operator 

^ + Mt) ■ C°°{X x S\£) — > C°°(X x S\£) (10) 

2 ' Recall that a linear first-order differential operator (|6]) is called a Dirac type operator if its principal symbol 
cta(x,£) satisfies the condition (cta{x, £)) 2 = ^2 ff 7 '* ' { x )(.j£kl where / is the identity operator in the fiber E x 
and the g-* k (x) are the (contravariant) components of the metric tensor (see [28]). The formal self-adjointness 
of A is understood in the standard sense as the condition that the identity (u, Av) = (Au, v) holds for any 
sections u, v € Cq°(X \ dX, E), where (•, •) is the inner product on L 2 (X, E) = L 2 (X, E,dvoV). Here dvol is the 
Riemannian volume element on X. 

3 ' Apparently, one derivative would suffice, but let us think big. 

4 -'We assume the circle S 1 to be obtained from the interval [0, 1] by gluing together the endpoints. 
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Figure 3: The manifold X x S\ obtained by gluing together the faces t = and t = 1 of the 
product X x [0, 1], and the operator + A{t). 

is elliptic, and the boundary conditions 

irc(u\YxS x ) = 0' where 7T£: fyxS 1 — ^ ^YxS 1 /^- ^ s the natural projection, (11) 

defined by the subbundle C, are elliptic for the operator (|l(jp , TTte closure (Jj + ^4(t)) £ o/ the 
operator ()10p /rom f/ie domain specified by conditions (jlip is an unbounded Fredholm operator on 
L?{X x S" 1 ^) with domain T>c consisting of the sections u 6 H 1 (X x S ,£) satisfying condition 

dm 



Now we are in a position to state the main theorem of the present paper. 
Theorem 2. One /ias 

B f{A(i) i(t) }=ind(^ + A(i) 



(12) 



The right-hand side of (fT2]) is the analytic index of the operator + A(t)) £ , i.e., the 
difference of dimensions of its kernel and cokernel, which can be expressed in topological terms 
by the Atiyah-Bott formula [H] (see also [22 s Sec. 20.3]). 

4 Proof of the Main Assertions 

Proof of Proposition [TJ Consider the operator 




d_ 

hi 



+ A{t) a 'V (t) ) ■■ C°°(X x S\£ ®£) ^ C°°(X x S\£ ©£). 
This is a total formally self-adjoint Dirac type operator on X x S 1 with symbol 
„ + t t \ - ( i£oI + o- A ( t) (x,£) 



(13) 



(14) 
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where £o is the momentum variable conjugate to t G S , and the operator ^ + A(t) is its 
chiral part. The subbundle £ = £®£c(£® £) YxS i satisfies a condition of the form © 
with respect to a% and hence specifies self-adjoint elliptic boundary conditions for 21. Indeed, 
the conormal vector to the boundary of X x S l at an arbitrary point (x, t) G Y x S 1 has the 
form n(x,i) = (0, n(x)), where n(cc) is the conormal vector to the boundary of X itself and 
£(x,t) = L (t)x) hence, for any 

V = t (v±,V2) G £(E,t)i 1" = t {wi,W2) G £(a;,t)> i.e., UljW2 3 UJ l>' u, 2 € L(t)x) 

we have 

(u,cra(x,t,n(x,t))io) = (u x ,(T A ( t )(a;,n(a;))ti;2) + (v 2 ,cr A ^(x,n(x))wi) = 0, 

because condition ([7]) is satisfied for A(i) and the bundle L(t). This, again by virtue of the results 
in [29] and [2S\ Chaps. 18 and 19], implies the claim of Proposition [T] first for the operator 21 
and then, as a consequence, for its chiral part + A(t). □ 

Proof of Theorem^ a. Without loss of generality, we assume that ^ Spec(A(0)£(o))- (Oth- 
erwise, one can replace the operator A(t) by A(t) + e with small real e, which changes neither 
the left- nor the right-hand side of (|12f) .) 

b. Also without loss of generality, we assume throughout the following that the subbundle 
L(t) is independent of the parameter t, L(t) = L(0) = L, t G [0, 1]. Indeed, let V(t) : Ey — > Ey 
be a family of unitary automorphisms of Ey such that V(0) = / and L(t) = V(t)L, t G [0, 1]. 
(One can readily construct such a family by solving the Cauchy problem V = [P, P]V, V(0) = /, 
where P = P(t) is the projection onto L(t) in Ey.) This family can be continued (by a homotopy 
to the identity mapping along the variable normal to the boundary) to a family of unitary 
automorphisms W(t): E — ^ E such that W(i)| y = V{t). Set 

A(t) = W-\t)A(t)W{t); 

then, obviously, 

TT _1 (l)ETA(0)(W'~ :l (l)i7)~ 1 = W^i^UA^U^Wil) = W -1 (1)A(1)W(1) = A(l), 
W~ l {l)UL = W~\\)UL{fS) = VF~ 1 (1)L(1) = L(0) = L; 

i.e., conditions of the form ([9]) are satisfied for the family {A(i)} and the constant family of 
subspaces L(t) = L if one takes the automorphism U = W~~ l {l)U . Furthermore, 

sf{A(t) m } = si{A(t) L }, (15) 

because the operators A(t) and A(t) are similar. Next, the family W(t) generates a bundle 
isomorphism W: £ — > £ , where the bundle £ over X x S 1 , by analogy with £, is obtained from 
the pullback of E to X x [0, 1] by clutching with automorphism U. The operator 

has the same index as ( + A(i)) £ and acts on the space of sections of £ satisfying the boundary 
condition associated with the subbundle C = W _1 | yxS1 £, for which C t = V(t)~ l L(t) = L for 
all t G S . Finally, the homotopy 

+ A(t) - \W- l (t) — {t)} A G [0, 1], (17) 
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in the class of Fredholm operators reduces the operator (fT6|) for A = to (Jj + A(t)) ^, so that 



ind(^ + A(t) 



i„d( ^ + A(t) 



which, together with (fT5|) . completes reduction to the case of a bundle L(t) = L independent of 
t. We omit the tilde over letters in what follows. 

c. In the proof, we need a family of operators on the infinite cylinder IxR. Let us describe 
it. The pullbacks of the bundle E from X to X x R and the bundle L from Y to Y x R will 
be denoted by the same letters E and L, respectively; this shall not lead to confusion. The 
coordinate on the line R will be denoted by t. For a, ft G R, we introduce the weighted spaces 
L 2 a JX x R, E) and H^JX x R, E) of sections u of E with finite norm 



u 



\u\ 



0,a$ 



L 2 (X,E) 

du(t) 2 



e 2at dt + 



W(t)\\h(x,E) e Wt dt 



1/2 



and 



dt 



+ 





0u(t) 




0* 



L 2 (X,E) 
2 



+ « 



H X (X,E) 



e 2at dt 



L 2 {X,E) 



+ u 



2 \ p m Hf 



1/2 



respectively. In particular, H^AX xl,£) C L? a ^(X xR,£). By 2) aj g we denoted the closed sub- 
space of H^,g (X x R, E?) consisting of the sections satisfying the boundary conditions determined 
by L; i.e., 

= {u£ Hl p (X xl,E): vr L -u = 0}. 

Let < 9 < 1. Se#i] 



(18) 



r(t,9) 



t < 0, 
< t < 
9 < t. 



For 7 G R, let 



L 2 Ql (XxR,E) 



(19) 



be the operator with domain Dq 7 (see Fig. Let us state a number of properties of the 
operators A(9, 7). 

Lemma 1. The operator A(9,j) is Fredholm for 9 such that 7 ^ Spec(A(9) l) , and ind A(9,'j) 
is a locally constant function of 9 on the set of such values of 9. 

Lemma 2. 1/7,7 ^ Spec(A(9)i) and'j > 7, then the difference ind A(9 ,^y) — ind A(9 , j) is equal 
to the number of eigenvalues (counting multiplicities) of the operator A(9)l on the interval (7,7). 

Lemma 3. indA(0,0) = 0. 

Lemma 4. ind A(l, 0) = ind(^ + 

The proof of Lemmas [THU will be given below. Now let us show that these lemmas imply 
the claim of the theorem. The spectral flow of the family {^(i)^} is given by Definition [T] for 
some partition = to < *i < ■ ■ ■ < tn+i = 1 and real numbers 71, . . . ,7n+i = 7l> and we can 
assume that 71 = (because we have assumed that ^ Spec( J 4(0)i)). Let rrij be the number of 



5 'Here, just as above and below, we for brevity omit the standard smoothing procedure eliminating the jumps 
of the derivatives (in the present case, for t = and t — 6) when describing the homotopies. 
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I + Mf>) 






XxR 



1 




r(t,l) 


e 




T (t,e) 










o e 




t 



Figure 4: The operator A(9,'y) = Mr + A(r(t, #)) on the infinite cylinder XxR. 



eigenvalues (counting multiplicities) of the operator A(tj)L in the interval between jj and 7j+i- 
It follows from Lemma [1] that the operator A(9, jj) is Fredholm for 9 G [tj_i,tj] and 



ind^4(tj,7j) — md^4(tj_i, 7j) = 0, j = 1,. . . ,n + 1. 



By Lemma [21 



ind^4(tj,7j + i) - ind.4(t,-, 7,) = raj sign^ - 7j+i), i = 1, 



, n. 



(20) 



(21) 



By summing relations (J20J) and f|21|) over all corresponding j, by adding the results, and by 
taking into account Lemma [3] and the relation 71 = 7 n = 0, we obtain 



ind.4(l,0) = ind.4(l,0) -ind„4(0,0) = ^m i sign( 7i - = sf{,4(i) L }. (22) 

3=1 



It remains to use Lemma HI The proof of Theorem [2] is complete. 
Now let us prove Lemmas HHH 



□ 



Proof of Lemma [TJ To prove that the operator A(9, 7) is Fredholm, it suffices to construct a 
regularizer, i.e., an operator 

TZ: Lg 7 (X x R, E) — > D 07 

such that the operators / — A(9, ^y)TZ and I — 1ZA(9, 7) are compact in the spaces Lq^(X x R, E) 
and 2)07, respectively. This can be done by the frozen-coefficients technique, standard in elliptic 
theory. (In our case, we "freeze" the variable t). To this end, for given r E [0, 1] and v € R, 
consider the operator 



dt 



+ A(t) : L 2 UU (X x R, E) — > L 2 UU (X x R, E) with domain ©„„. 



(23) 



The operator (f23j) is invertible provided that v ^ Spec( J 4(r)z / ). The inverse operator R u (t) is 
given by the formula 

[R v {r)u]{t) = f e^iip + A(r) L )- x u{p) dp, u G L 2 W [X X R, £), (24) 

V 27T Jl mp = u 



9 



where u(p), Imp = u, is the Fourier transform of u with respect to the variable t. Consider 
a finite cover {Uj}j =1 of [0,1] by open intervals such that Uj ^ Spec(A(r(t, 9))l) for t G Uj 
for some real numbers Vj\ let Uq = (— oo,0), uq = 0, U s+ \ = (l,oo), and v s +i = 7, and let 
1 = Ylj=o ipj be a smooth partition of unity subordinate to the cover of the line R by the sets 
Uj, j = 0, . . . , s + 1. Then 1Z can be defined by the formula 



+i 

r(t,ey 



[nu\{t) = Y,Ut)[[n V] {T){^u)\{t)} 

j=0 



(Note that 1Z is well defined as an operator from L^{X x R, E) to So 7 , because the operator 
of multiplication by ipj is continuous from Lq^(X x R, E 1 ) to 1% v (X x R, E) and from S)^^ to 
2) 07- For j / 0, s + 1, this follows from the compactness of the support of ipj; for j = 0, from 
the fact that uq = and ipo(t) = for i > 0; for j = s + 1, from the fact that za^+i = 7 and 
-(/>,j+i(i) = for t < 0.) Now a straightforward computation shows that 



^(0,7)tt«(t)=«(f) + 0[(^(f)^ 

j=o LLV 



(*) 



Since the functions ~^-{t) and ^(t,9) are compactly supported, it follows from standard facts 
about embeddings of Sobolev spaces that the second term on the right-hand side defines a 
compact operator on Lq^(X xR,i?), In a similar way, one can study the product 1ZA(6,j). 

The local constancy of the index of A(9,j) as a function of 6 follows from the fact that 
this operator continuously depends on 9 in the operator norm as an operator from £>07 to 
Lq 7 (X x R, E). The proof of Lemma [1] is complete. □ 

Proof of Lemma [3j This lemma is the special case of the invertibility of the operator (|23p for 
v = and r = 0. □ 

The proof of Lemmas[2]and[3]is based on the localization method (the index locality principle; 
see [26} Theorem 4.10] and also [23-25] and the references therein). Having in mind our goals, 
let us state the claim of Theorem 4.10 in [26J for the simplest special case. 

Let N\,N2 C M be disjoint closed subsets of a manifold M, and let D: H\ — > %2 be a 
bounded Fredholm operator acting on some Hilbert spaces of sections of bundles over M. Next, 
let x: M — > [0,1] be a smooth mapping such that N\ C / _1 (0) and N2 C /~ 1 (1), and let 
C C C°°(M) be a subalgebra consisting of functions constant on N\ and on N2 and containing 
all functions of the form ip(x) = (p(x(x)), where (p is a smooth function on [0, 1]. Suppose that 
the commutator of D with the operator of multiplication by any function in C is compact. Then 
the index increments arising from changes of D on N% and N2 preserving the Fredholm property 
and the compactness of commutator i^l are independent: 



where 



A ^ is the index increment occurring if the operator is changed only on N\ . 
A jv 2 is the index increment occurring if the operator is changed only on N2 ■ 



6 'The changes may affect not only the operator itself but also the spaces on which it acts and even the very 
manifold (e.g., cutting away some parts and pasting another ones); all these changes should occur strictly inside 
the corresponding set Nj, and everything on M \ Nj should remain unchanged. 
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• A j /v 1 uAf 2 is the index increment occurring if the operator is simultaneously changed both on 
Ni and N 2 . 

The practical application of the localization method to the proof of Lemmas [2] and H] imple- 
ments the following idea. We wish to compute how the index of some operator D changes for 
a given change of the operator on a set N\, but it is difficult to compute the index increment 
owing to the complicated structure of D outside N\ . Let us modify the operator D on some set 
N2 disjoint with N\ so as to obtain an operator D of simpler structure whose index increment 
under the given change on JVi can be computed. This increment coincides with the desired 
increment for the original operator. 

Proof of Lemma EJ Take X x R for the manifold M, the set {t > 2} for Ni C M, the set {t < 1} 
for N2 C M, and the algebra of infinitely differentiable functions (p(t) of t € R constant on N\ 
and on iVjj for the function algebra C. The original operator D is the operator A(9, 7), which 
we treat as a bounded Fredholm operator on the spaces 

D = A(fl l7 ):J)» y -f4 y (XxR I £) I (25) 

and we need to compute the index increment for this operator if 7 is replaced by 7. Note 
that the commutator of the operator (|25p with a smooth function ip £ C is the operator of 
multiplication by the compactly supported function tp'(t), which is compact as an operator from 
Do-y to Lq^(X x M 1 ,^), so that we are just in a position to use the localization method. The 
replacement of 7 by 7 changes the operator D only on the set N\. (The expression specifying 
the operator and the boundary conditions remain the same, but the spaces where the operator 
acts are changed, the change being solely concerned with the admissible growth of functions as 
t — > +00; i.e., in particular, the restriction of these spaces to M \ Ni is unchanged at all.) Now 
let us replace D by the operator 

D = ^ + A(9):T)^^ L 2 0j (X xR,E). (26) 

This operator differs from D (both in the differential expression and in the spaces where it acts) 
only on N 2 - Thus, it suffices to compute the index increment for this operator under the change 
on Ni the same as for D. The operator D is invertible, so that mdD = 0. The change of this 
operator on N\ results in the operator 

^ + A(«):%->I^xR,£); (27) 

thus, it remains to compute the index of the latter. For this computation, it is convenient to 
treat the operator (|27p as an unbounded Fredholm operator on L^~(X x R, E) with domain J)^. 
Then the adjoint operator has the form — ^ + A{9) and acts on the dual space _^(X x R, E) 
with domain 2)_ 7 _^. The elements of the null space of the operator (|27p should have the form 
v(x)e~ >yt , where A is an eigenvalue of A(9)l and v(x) is one of the corresponding eigenfunctions. 
The condition that these elements belong to the weighted space L?~(X x R, E) implies that 
7 < A < 7. Thus, the dimension of the null space is equal to the number of eigenvalues (with 
regard of multiplicity) of the operator A(9)l in the interval (7,7). The elements of the null 
space of the adjoint operator should have the form v(x)e xt and belong to the weighted space 
£?. 7 _ 7 (-X" x R> E)). It follows that 7 < A < 7, but this is impossible, because 7 > 7. Thus, the 
null space of the adjoint operator is trivial, the index of the operator (|27|) coincides with the 
dimension of its null space, and we arrive at the assertion of Lemma [2j □ 
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Proof of Lemma [^} We need to prove that the operators 

.4(1,0): Doo — >L 2 (XxR,E) (28) 

and 

^ + i(t)J^ £ ^L 2 (IxS^) (29) 

have the same index. We assume (this can always be achieved by a homotopy) that A(t) = 4(1) 
for t > 1 — 2e and A(t) = A(0) for t < 2e for some given e > 0. Set 

Ni = {te (-oo,e] U [1 -e,oo)}, N 2 = {t € [2e, 1 - 2e]}. 

For the function algebra C we again take the algebra of infinitely differentiable functions (p(t) 
of t constant on N\ and N 2 . The operator f)29[> can be obtained from the operator (|28p by 
the following change on N\: one cuts away and disposes of the half-cylinders X x (— oo,0) 
and X x (l,oo), and the faces X x {0} and X x {1} of the remaining product X x [0,1] 
are glued together, the bundle E giving rise to the bundle £ via the clutching automorphism 
U : E\t=o — > E\t=i. 

We should show that the index increment for this change of the operator on N% is zero. To 
this end, we replace the operator ([28]) by the operator 



4(0,0): S) 00 — > L 2 (X x R,E). (30) 

The operator (|30p will differ from the operator (|28p only on N 2 if we rewrite the former in the 
equivalent form 

M00)-l^ + m ***** (31) 



where it is assumed that the bundle in which the operator (|3ip acts is obtained from E\ t< i and 
E\t>i by the standard clutching construction at t = ^ with the automorphism 



U : E\ t=1 / 2 +o — >■ E\t=i/2-o- 

(The passage from (|30p to pip is essentially none other than rewriting the operator 4(0, 0) for 
t > § i n "new coordinates" in the fibers of -E 1 .) 



Now let us change the operator pop written in the form pip on JVi in the same way as we 
have earlier changed the operator (|28p . The resulting operator on X x Si has the form pip , and 
the bundle in which it acts is obtained from E by clutching construction with the automorphism 

U:E\ t=0 — > E\ t= i, U ■ E\ t=1 / 2+0 — > E\t=i/2-o- 

It is easily seen that this bundle is isomorphic to the pullback of E on X X Si, and the resulting 
operator itself is none other than (Jj + 4(0)) L ; its index is zero, because it is invariant with 
respect to rotations along S . The index of the operator ()28p is zero as well (it is invertible), so 
that the index increment is zero, and the proof of the lemma is complete. □ 

Proof of Theorem{l[ Theorem [2] shows that the spectral flow of the family ([2]) obeys the localiza- 
tion principle: any modifications applied to the Dirac operator (JTJ) in the planar domain X and 
to the function B occurring in the boundary conditions automatically lift to X x Si becoming 
modifications of + D±\ the latter enjoy the index locality principle [26|, Theorem 4.10], and 
the index is equal to the spectral flow by Theorem [2j 

The localization principle permits one to split the domain with holes into parts with fewer 
holes. Let us show this by example. Figure [SJ left shows a domain X with two holes. Let us 
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Figure 5: Cutting a domain with holes into pieces. 




Figure 6: Cutting a domain without holes into pieces. 



reduce the computation of the spectral flow of a family of Dirac operators with local boundary 
conditions in X to the corresponding computation for domains with one hole. Let N± be the 
set dashed in Fig. [5J left, and let N2 be the complement to a small neighborhood of N\. Let us 
change the domain inside N\ as shown in Fig. [SJ right, so that the original domain becomes two 
domains with smooth boundary. The function B{x) occurring in the boundary conditions can 
be extended by continuity as a nonvanishing real-valued function to the newly arising boundary 
arcs inside N\, because the sign of B(x) is the same on the entire outer boundary of the original 
domain. Thus, the domain splits into two unrelated parts, and to prove that the spectral flow 
of the family of Dirac operators in the original domain is equal to the sum of spectral flows 
corresponding to the two new domains, one should show that the increment of the spectral flow 
under this modification of the domain is zero. To this end, we use the localization method. Let 
us change the original family by changing the domain in N2 (so that the resulting domain has the 
form shown in Fig. [UJ left) and by extending B{x) by continuity as a nonvanishing function to 
the newly arising boundary arcs inside N2. The spectral flow of the new family is zero before as 
well as after the modification shown in Fig. because the domains in this figure are contractible 
and the gauge transformation \i in these domains is homotopic to the identity transformation. 
Thus, cutting the domain into pieces reduces the problem to the case of domains with one hole, 
for which formula ([3|) was proved in [20J. (Needless to say, one can prove it directly with the 
use of Theorem [21 but there is no need to do this, and we omit the corresponding computations 
for lack of space.) □ 

5. CONCLUSIONS 

Let us discuss possible physical realizations of nonzero spectral flow. We start with the case 
of graphene. One has to keep in mind that there are two Dirac electron subsystems in graphene 
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(two valleys) and, generally speaking, scattering at the edges mixes the two valleys [HE]- This 
is not the case, however, if an energy gap in the electron energy spectrum opens smoothly when 
reaching the edge. At the chemical functionalization of the edges this is, indeed, the case, since 
electronic structure is modified in a sufficiently large region of space [32]. As a result, intervalley 
scattering is negligible and we have the boundary condition ([T]) suggested first by Berry and 
Mondragon [19]. A detailed microscopic derivation of the boundary condition starting with a 
discrete lattice model has been done in Ref. [31] (see also Chapter 5] and references therein). 
The sign of the constant B is determined by the sign of the mass term in the Dirac equation 
and is dependent on the distribution of chemical groups along the edge. One can hope that if we 
prepare graphene rings, then in some specimens the signs of B will be opposite at the internal 
and external edges of the ring, which is necessary for nonzero spectral flow. However, it is hard 
to reach in a controllable way. 

Probably, topological insulators are more promising in this sense. First, two-dimensional 
massless Dirac fermions are realized at the surface of three- dimensional topological insulators, 
such as Bi2 Se3, only one Dirac cone arising [8,9]. To open the gap, one has to cover the surface 
by a magnetic layer, the sign of the gap being determined by the direction of magnetization [8,9]. 
This opens a way to manipulate the sign of the constant B. Second, two-dimensional massless 
Dirac fermions can be realized in a layer of HgTe confined between two layers of CdTe, at a 
certain critical thickness of the layer [8l[9] . Recently, such an opportunity has been demonstrated 
experimentally [33]. If the thickness of the layer varies smoothly in space oscillating near the 
critical value, one can reach both positive and negative values of B. Currently, this opportunity 
to create nonzero spectral flow looks the most promising. 
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